Doping dependence of thermopower and thermoelectricity in strongly correlated 

materials 
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The search for semiconductors with high thermoelectric figure of merit Q, Q, Q has been greatly 
aided by theoretical modeling of electron and phonon transport, both in bulk materials Q and 
in nanocomposites |^ 0|. Recent experiments 0, |E have studied thermoelectric transport in 
"strongly correlated" materials derived by doping Mott insulators, whose insulating behavior with- 
out doping results from electron-electron repulsion, rather than from band structure as in semi- 
conductors. Here a unified theory of electrical and thermal transport in the atomic and "Heikes" 
limit is applied to understand recent transport experiments on sodium cobaltate and other doped 
Mott insulators at room temperature and above. For optimal electron filling, a broad class of 
narrow- bandwidth correlated materials are shown to have power factors (the electronic portion of 
the thermoelectric figure of merit) as high at and above room temperature as in the best semicon- 
ductors. 

PACS numbers: 74.25.Fy, 73.50.Lw, 71.10.Fd 



The efficiency of a thermoelectric material for refriger- 
ation or electrical generation is determined by its ther- 
moelectric figure of merit 



ZT 



(1) 



where T is the temperature, a is the electrical conductiv- 
ity, S the thermopower and k the thermal conductivity. 
The same transport coefficients that determine a mate- 
rial's utility for thermoelectric applications are also of 
importance for fundamental materials science; for exam- 
ple, thermopower measurements are frequently used to 
measure the sign of the dominant charge carriers in a 
material. Most thermoelectrics currently in use Q are 
bulk semiconductors such as Bi2Te3, with ZT « 1 at 
room temperature. 

In recent years, strongly correlated oxides have been 
studied as candidate thermoelectric materials. Research 
in this direction was sparked off by the observation 
of high thermopower (~ 125iiV/ K) at room tempera- 
ture coupled with a low electrical resistivity in the ma- 
terial Na2;Co02 d, IE 1^' Other materials including 
SrT, Lai -irTiO.s and several superconducting cuprates 

were found to have promising values of thermopower 
at room temperature. These materials are strongly cor- 
related with an interaction energy scale U that is much 
larger than the bandwidth set by the hopping t. Theoret- 
ical calculations on such materials are difficult and often 
specific only to certain ranges of parameters 0, 0| . 

The "power factor" TS^a appearing in the numerator 
of ZT is determined by properties of the charge carriers 
(electrons or holes) , while k is dominated by phonons for 
most viable thermoelectric materials. Although applica- 
tions of oxides as thermoelectric materials are currently 
limited by their relatively high thermal conductivity, re- 



cent advances in reducing phonon thermal conductivity 
in semiconductor thermoelectrics 0, 0, Q| should lead to 
more dramatic gains in oxides, where the original ther- 
mal conductivity is higher. Hence, understanding which 
oxides can have power factors as high or higher than the 
best semiconductors is an important step toward realiz- 
ing the potential of this material class. 

In this paper, we extend the calculation of the ther- 
mopower of correlated systems in the atomic limit 15, 16| 
to the electrical and thermal conductivities and obtain 
the power factor and electronic ZT, as functions of dop- 
ing and temperature. Although the optimum band struc- 
ture for thermoelectric efficiency in a semiconductor was 
determined by Mahan and Sofo Q, the dominant role of 
electron-electron interactions in correlated materials re- 
quires entirely different theoretical methods. The power 
factor is found to display maxima at doping values of 5% 
and 88 % from the Mott insulator consistent with recent 
experiments on Naa;Co02 [3- We discuss the effect of 
multiple orbitals in closing and comment on the search 
for high-efficiency thermoelectric oxides. 

We describe a strongly correlated system by the Hub- 
bard model, which consists of a lattice with electrons 
capable of hopping from a site, with a single orbital, to 
its nearest neighbors with an energy (hopping parame- 
ter) t. Double ocuupancy incurs an energy cost of U per 
site. This surprisingly complex model has been used to 
study the physics of superconductivity, spin and charge 
ordering, and the Mott transition seen in a variety of 
oxides. The Mott insulator appears at the occupancy of 
one electron per site. At temperature T, the atomic limit 
is defined as the approximation t ^ {kBT,U). Systems 
with narrow bandwidths (e.g., Naa;Co02) are best suited 
for this type of analysis. In this limit, cr, S and k are 
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given by the following expressions to leading order 



^/3(C/-3m) 



S 
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/3/i 



(2) 



(3) 



and 



(/3C/)2e-''(^-3'') 



(4) 

Here /? = l/ksT and ^ is a doping and temperature 
independent, but material dependent, factor given by 



A = 



Vh? ' 



(5) 



where a is the lattice spacing, y, the volume of the unit 
cell and r, the transport relaxation time. 7 is a geometry 
dependent constant equal to 1 for a square lattice and 2 
for a triangular lattice. 



= ( 1 + 26^3^ + e-/3(c/-2M) 



(6) 



The chemical potential /i is expressed in terms of the 
doping p as 



p-l + ^J{p-lf+p{2-p)e 
(2 - p)e-/3'^ 



-/3(7 



(7) 



p = 1 is the Mott insulator and p — Q and p = 2 cor- 
respond to completely empty and full bands respectively 
and are both band insulating limits. The second term in 
the expression for the thermal conductivity k is of sub- 
leading order at high temperature {ksT U) but not in 
the Heikes limit ksT <C U, where it is the leading term. 
It should be noted that both a and k depend on a re- 
laxation time T as in Drude theory, while S and ZT do 
not. There are further simplifications in the Heikes limit 
ksT-^U: 



S 



■log 



2(1 -P) 



P(l-P), 



(8) 



(9) 



The above expressions are for p < 1 and can be ex- 
tended to p > 1 by making the substitution p — > 2— p. We 
note that the thermopower has no temperature depen- 
dence and is independent of material parameters in this 
limit; it changes sign by going through a resonance at the 
Mott insulator p = 1- Experimental data on Naa;Co02 is 
compared with theory in Fig. 1. The electronic ZT is di- 
vergent close to the Mott insulator but adding a phonon 
contribution to k keeps it finite. This is illustrated with 
the material parameters for Naa;Co02 as an example 
in Fig. 2 




FIG. 1: The experimentally measured thermopower (aster- 
isks) at different values of doping for Naa;Co02 at 300 K. The 
theoretical Heikes curve is provided alongside for comparison. 
It can be seen that the experimental values are close to the 
theoretical ones, showing that the Heikes limit is a good ap- 
proximation at 300 K. The values are systematically lower 
since the Heikes limit is approached from below in tempera- 
ture. 
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FIG. 2: at 300 K for different values of doping for 

Naa;Co02 (with t — llOK) assuming different values of the 
the phonon thermal conductivity Kph (assumed to be doping 
independent). Kph adds to k from Eq. 10. 

The relaxation time t ~ 2 x 10^^^ s for Naj-CoOo is 
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extracted from the conductivity measurements |9| . a; = 
corresponds to the Mott insulator p = 1- A plot of the 
theoretical TS'^a for Naa;Co02 is shown in Fig. 3. The 
value of TS'^a (« 1 — 2Wm~^K~^ is indeed observed in 
experiments. The main result of the above calculation 
is a robust maximum in the power factor at dopings of 
p — 0.12 and p = 1.88 and secondary maxima close to 
the Mott insulator at p = 0.95 and p = 1.05. These 
maxima arise from competition between the (large) ther- 
mopower and (small) conductivity near the Mott {p = 1) 
and band {p = and p = 2) insulators, p = 1.88 cor- 
responds to x=0.88 for Naa;Co02, where a peak in the 
power factor has indeed been observed in recent experi- 
ments over a range of temperatures . Although we have 
used Naa;Co02 to enable comparison with experimental 
data, the fact that the main features are reproduced in 
the Hubbard model shows that any material with roughly 
similar values of t, U, and r will have the same features. 



materials, such as La^jSri-a-TiOa (for which polaronic 
effects are known to be important), close to the band 
insulator [loj . 

We finally examine the effects of multiple orbitals 
per site on thermoelectric transport. This model is 
applicable to certain strongly correlated materials like 
chromium and magnesium based oxides We con- 

sider a multiple-orbital Hubbard model with N degen- 
erate orbitals per site and a single value of the hopping 
parameter (for concreteness), and extend the calculation 
for the thermopower [3| to include a, k, ZT and TS^a. 
There is a correlation energy of U for each pair of carriers 
occupying the same site. In the Heikes limit, we obtain 
the following expressions for the transport coefficients: 

(T=— ^(1-p), (12) 



^Theoretical estimate for Na^CoO^ at 300 K 
x^1-p 
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FIG. 3: The power factor (TS'^a) as a function of doping 
for Naa;Co02 as predicted by theory. There are maxima at 
p « 0.12, p « 0.95, p « 1.05 and p « 1.88. The last two 
values can be realized in the material in principle at doping 
values X = 0.05 and x = 0.88. The positions of the peak are 
independent of the material parameters in this calculation 
and a peak at a; = 0.88 has indeed been seen in experiments 
on Naa;Co02. 

The location of the doping peak depends only weakly 
on temperature and for Na2;Co02 and shifts appreciably 
only at r > 700 K. The theoretical position of the peak is 
independent of material dependent parameters. In reality 
however, factors like the effect of phonons on electron 
transport (phonon drag and polarons), orbital degrees 
of freedom, etc. could change the location of the peak. 
As mentioned earlier, the above theory initially assumed 
the atomic limit and hence is best applied to narrow- 
bandwidth systems. For other correlated systems, one 
might need to consider terms beyond leading order in 
the calculation of transport coefficients. However, the 
existence of a peak in the power factor is quite general 
and has been observed in other candidate thermoelectric 



(13) 



(14) 



TS'^a for an arbitrary system at a fixed temperature 
T is shown in Fig. 4. There are still strong maxima 
as a function of doping, though displaced slightly from 
the single-orbital [N — 1) value. This is due mainly 
to the logarithmic dependence of S on N . The overall 
magnitude increases for large N ^ as the presence of more 
hopping channels increases the conductivity. 




FIG. 4: The value of the power factor at a particular tem- 
perature for an arbitray system, as a function of doping for a 
multiple-orbital Hubbard model. All material dependent pa- 
rameters are assumed to be the same for the different values 
of A'^. Larger A'^ implies more conducting channels giving a 
larger conductivity and hence a larger value of TS^a. The 
peaks shift with increasing A'^ but not by much owing to a 
weak logarithmic dependence of the thermopower on A'^. 
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To conclude, we have investigated the thermoelectric 
properties of strongly correlated systems by studying the 
Hubbard model. We have shown that in the absence 
of phonon heat conduction, ZT can diverge at values 
of doping close to the Mott and band insulators. The 
calculated power factor and thermopower agree with ex- 
periments on Na2;Co02. We have in particular demon- 
strated the existence of a doping maximum 88 % away 
from the Mott insulator, also been seen in experiments 
and have investigated the effect of multiple orbitals on it. 
The prediction of a large power factor and the presence 
of a doping maximum should aid in searches for better 
thermoelectric oxides. 
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from DOE. 



[1] H. J. Goldsmid, Electronic Refrigeration (Pion, London, 
1986). 

[2] D. M. Rowe and C. M. Bhandari, Modem Thermo- 
electrics (Reston Publishing, 1983). 



[3] F. J. DiSalvo, Science 285, 703 (1999). 
[4] G. D. Mahan and J. O. Sofo, Proc. Natl. Acad. Sci. 93, 
7436 (1996). 

[5] M. S. Dresselhaus et al.. Mater. Res. See. Symp. Proc. 

478, 55 (1997). 
[6] G. Chen, Phys. Rev. B 57, 14958 (1998). 
[7] I. Terasaki, Y. Sasago, and K. Uchinokura, Phys. Rev. B 

56, R12685 (1997). 
[8] Y. Wang, N. S. Rogado, R. J. Cava, and N. P. Ong, 

Nature 423, 425 (2003). 
[9] M. Lee et al.. Nature Materials 5, 537 (2006). 
[10] T. Okuda, K. Nakanishi, S. Miyasaka, and Y. Tokura, 

Phys. Rev. B 63, 113104 (2001). 
[11] S. D. Obertelli, J. R. Cooper, and J. L. Tallon, Phys. 

Rev. B 46, 14928 (1992). 
[12] A. Oguri and S. Maekawa, Phys. Rev. B 41, 6977 (1990). 
[13] B. S. Shastry Phys. Rev. B 73, 085117 (2006). 
[14] R. Venkatasubramanian, E. Siivola, T. Colpitts, and 

B. O'Quinn, Nature 413, 597 (2001). 
[15] G. Bcni, Phys. Rev. B 10, 2186 (1974). 
[16] S. Mukerjee, Phys. Rev. B 72, 195109 (2005). 
[17] B. Kumar and B. S. Shastry, Phys. Rev. B 68, 104508 
(2003). 

[18] W. Kobayashi, L Terasaki, M. Mikami, and R. Puna- 
hashi, J. Phys. Soc. Jpn. 73, 523 (2004). 



